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Abstract.  Study  of  the  deviations  of  P  and  S 


travel  times  from  the  J-B  tables  at  teleseismic  distances 
has  shown  that  there  are  regional  differences  in  travel 
time.  Both  P  and  S  are  early  in  the  central  and  eastern 
United  States,  late  in  the  western  United  States.  The 
differences  have  a  range  of  about  three  seconds  for  P 
and  eight  seconds  for  S. 

It  can  be  deduced  from  the  relation  between  the 
travel  time  residuals  (1)  that  the  change  in  shear  velocity 
is  approximately  one  and  one-quarter  times  the  change  in 
P  velocity,  (2)  that  the  observations  imply  a  difference 
in  Poisson's  Latio  between  the  two  regions,  and  (3)  that 
a  model  in  which  the  shear  modulus,  u,  alone  varies,  the 
compressibility,  k,  remaining  sensibly  constant,  fits  the 
data  best.  It  can  be  shown  also  that  the  differences 
between  the  P  travel  time  residuals  and  the  gravity 
anomalies  in  the  central  and  western  United  States  are 
not  consistent  with  the  Birch  relation  between  velocity 


•PMUPFISC  &mx  zi**» 


m 


» 


m 


,£F 


_  2  _ 


INTRODUCTION 

One  of  the  early  results  of  the  study  of  the 
deviations  of  observed  P  travel  times  from  the  Jeffreys- 
Bullen  table  times  by  Cleary  and  Hales  (1963,  1965,  1966) 
was  that  there  was  a  significant  component  of  the  deviations 
which  was  regionally  dependent.  Recently  Doyle  and  Hales 
(1967)  have  studied  the  deviations  of  S  travel  times  from 
the  tables  using  the  same  method  of  analysis  as  Cleary  and 
Hales.  The  S  station  residuals  for  North  America  ranged 
over  eight  seconds  in  contrast  to  the  P  residuals  for  which 
the  range  was  only  two  to  three  seconds.  The  P  and  S 
residuals  at  the  same  station  were  correlated,  the  correlation 
coefficient  being  0.75.  Figure  1  shows  a  plot  of  the  S 
residuals  against  the  P  residuals  after  adjustment  to  make 
the  S  residuals  zero  in  the  same  region  as  the  P  residuals 
(see  Doyle  and  Hales,  1967).  A  straight  line  fitted  to 
the  data  using  a  method  due  to  York  (1966), had  a  slope  of 
3.72  ±  0.43.  This  line  is  shown  in  Figure  1.  It  is  the 
purpose  of  this  paper  to  explore  the  implications  of  the 
rather  large  ratio  of  the  S  residuals  to  those  for  P. 

Residuals  of  this  order  of  magnitude  can  only  arise 
as  a  result  of  regional  differences  in  the  structure  of  the 
upper  mantle.  Let  us  compare  portions  of  two  upper  mantle 
structures,  on^  with  P  and  S  velocities,  aQ  and  0fl  respectively, 
and  the  other  with  velocities  afl  -  5a,  0O  -  60,  and  assume 
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that  the  paths  through  the  regions  of  difference  each  have 

length  D.  Then,  to  first  order  at  least,  the  difference  of 

the  P  travel  times  along  the  two  paths,  6t_,  will  be  , 

p  a. 2 

0  2 

whereas  6te  will  be  .  The  ratio,  6tc/6t_,  is  •  ^L_  . 

S  s  2  S  ?'  6a  o  2  • 

o  p  o 

In  deriving  the  above  results  the  effects  of  the  differences 
in  path  geometry  arising  from  possible  differences  in  the 
ratios  6a/ao ,  and  6 8/ S 0  have  been  neglected.  In  general  the 
effects  of  changes  in  the  path  geometry  will  be  small 
compared  with  the  effects  of  changes  in  the  velocities.  We 
have,  therefore, 


6tS  =  68  #  a0 

ZtZ  Ta  ’  a  : 


=  3.72. 


a 

Assuming  that  =  /3,  corresponding  to  a  Poisson's  ratio  of 
0.25,  we  find  68  =  1.24  6a. 

If  we  assume  that  a0  =  8.2,  8n  =  8.2,  /I  =  4.734,  it 
follows  that  a  =  4.38  +  0.818.  Figure  2  s^ows  a,  a/6  and  a, 
Poisson's  ratio,  as  functions  of  8  for  this  modelt  It  should 
be  noted  that  the  relation  between  a  and  6  given  above  applies 
only  to  those  portions  of  the  upper  mantle  in  which  the  travel 
time  differences  arise. 


MODEL  CALCULATIONS 

Many  studies  of  the  structure  of  the  upper  mantle 
have  been  carried  out  based  upon  the  dispersion  of  surface 

waves.  In  fact  general  acceptance  of  the  Gutenberg  low 


♦Footnote:  Similar  conclusions  to  those  of  this  paDer  follow 


for  other  values  of  8Q. 
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velocity  layer  in  the  upper  mantle  came  from  such  studies. 

Many  of  these  were  studies  of  average  structure  over  large 
areas,  while  others  dealt  specifically  with  local  areas. 

Of  these  the  two  most  relevant  to  the  present  discussion  are 
those  by  Alexander  (1963)  of  the  western  United  States  and 
McEvilly  (1964)  of  the  central  United  States. 

We  have  calculated  travel  times  at  distances  of  20  to 
80  degrees  for  a  number  of  upper  mantle  S  velocity  distributions 
based  on  surface  wave  and  other  studies.*  In  most  of  the  cases 


the  distributions  have  been  adjusted  to  fit  a  Gutenberg  S 
velocity  structure  at  some  depth  depending  on  the  range  over 
which  the  original  model  was  defined.  The  Gutenberg  velocities 
used  were  those  given  by  Press  in  the  Handbook  of  Physical 
Constants  (Press,  1966) . 

Table  1  lists  the  sources  of  the  models,  while  Table  2 
gives  the  velocity  distributions.  In  calculating  the  travel 


•  ,  V'  •  ^ 

times  the  Mohorovicic  law 


r  .’Ik 

,  v  =  v.  —  , 

i  r  ' 


log  v j  -  log  v, 

where  k  =  -» — - — - 1  , 

log  r  -  log  r2 


was  used  between  the  tabulated  points.  Figure  3  shows  the 
differences  in  travel  time  of  the  models  from  those  for  the 
Gutenberg  S  velocity  distribution  over  the  distance  range 
25°  to  85°.  It  will  be  noted  that  the  differences  decrease 
by  about  20%  from  30°  to  80°.  This  is  characteristic  of 

the  effects  of  P  and  S  velocity  changes  in  the  upper  mantle 

^Footnote:  The  models  were  selected  to  show  the  effects  of  a 
wide  range  of  upper  mantle  structures.  They  are  not 
comprehensive.  For  example,  travel  times  are  not  shown 
for  Miss  Lehmann's  models.  Teleseismic  travel  times  for 
these  models  lie  within  the  range  shown  in  the  figures. 
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and  indicates  that  station  anomalies  should  show  a  small 
decrease  with  distance.  For  the  McEvilly  Love  wave  model 
we  have  plotted  as  a  dashed  curve  the  diftererce  of  ScS 
travel  time  between  the  model  and  the  Gutenberg  distri¬ 
bution.  This  difference  changes  slowly  from  95°  to  25°. 

The  station  anomaly  for  ScS  is  always  less  than  for  S, 
ranging  from  94%  at  80°  to  77%  at  40°  of  the  S  station  anomaly 
at  the  same  distance. 

The  differences  in  travel  time  shown  in  Figure  1 


range  over  14  seconds  and  are  reasonably  consistent  with 
the  range  of  the  S  station  anomalies.  It  should  be  noted 
that  McEvilly  gave  two  S  velocity  distributions  for  the 
central  United  States,  one  based  on  Rayleigh  wave  dispersion, 
the  other  on  Love  wave  dispersion.  The  differences  between 
the  times  for  the  McEvilly  Love  wave  model  and  Alexander's 


western  United  States  model  average  10  to  12  seconds,  about 


3076  less  than  the  maximum  difference  of  8  seconds  in  the 
observed  S  times  (which  are  determined  for  one  end  of  the 
path  only) .  The  difference  between  the  times  for  the  McEvilly 
Rayleigh  wave  model  and  Alexander's  western  U.S.  model  are 
positive  and  about  one  second.  McEvilly  ascribed  the  differences 
between  the  Rayleigh  wave  and  Love  wave  models  to  anisotropy,  i.e. 
to  differences  in  velocity  of  SH  and  SV  waves.  For  the  travel 
time  anomalies  upon  which  this  discussion  is  based,  the 
upper  mantle  portion  of  the  path  was  in  most  cases  within 
30°  of  the  vertical,  i.e.  the  major  component  of  the 
particle  motion  was  in  the  horizontal  plane.  It  would  be 
expected  that  the  travel  time  anomalies  would  be  consistent 


- 


-J 


with  the  Love  wave  model  and  r.ot  the  Rayleigh  wave  model. 
Nevertheless,  the  discrepancies  shown  by  the  surface  wave 
dispersion  studies  in  the  United  States  are  puzzling  and 
merit  further  investigation. 

Figure  4  shows  that  changes  of  the  velocity  distribution 
in  the  upper  mantle  produce  effects  of  different  character 
from  those  resulting  from  changes  in  the  lower  mantle.  Plotted 
in  this  figure  are  the  differences  of  the  travel  times  between 
the  Gutenberg  model  and  a  number  of  other  models,  namely  the 
MuEvilly  Love  wave  model,  the  Jeffreys  model,  the  Anderson 
model  and  a  model  l  ased  on  the  Gutenberg  IV  model  of  MacDonald 
and  Ness  (1961)  (in  which  the  velocities  in  the  lower  mantle 
have  been  decreased) .  It  will  be  noted  that  the  differences 
for  all  models  other  than  the  Gutenberg  IV  model  vary  only 
one  or  two  seconds  over  the  teleseismic  distance  range, 
whereas  those  for  the  Gutenberg  IV  model  increase  systematically 
with  distance.  Comparison  of  this  figure  with  the  observed 
travel  times  of  Doyle  and  Hales  (1967)  and  Jeffreys  (1966) 
supports  the  conclusion  of  Landisman,  Sat6  and  Nate  (1965) 
that  changes  of  the  S  velocity  distribution  in  the  lower 
mantle  such  as  were  suggested  by  MacDonald  and  Ness  are  not 
permitted  by  the  observations. 

POISSON'S  RATIO 

It  can  be  shown  that  the  observed  relation  between 
the  P  and  S  residuals  is  not  consistent  with  a  constant 


Poisson's  ratio  in  those  parts  of  the  upper  mantle  for  which 

Ot  —  o  o  2 

the  regional  differences  occur.  If  a,  given  by  _L-* — .  , 

2(a*  -  3  ) 

were  constant,  then  a/6  would  be  constant,  6a/au  equal 
to  63/30*  and  the  rat.i.o  of  the  residuals  equal  to  a0/30, 
i.e.  1.7  to  1.8.  The  95%  confidence  limits  for  the  slope 
of  the  regression  line  of  S  and  P  residuals  are  2.86  and  4.58, 
and  thus  the  possibility  of  constant  a  is  remote. 

The  only  one  of  the  elastic  constants  which  can  be 
determined  independently  of  possible  changes  of  density  is  a. 
Expressed  in  terms  of  6a 

a  «  ^jl  -  . —  j  to  the  first  order, 

assuming  as  before  that  a0  ■  /580.  In  terms  of  5tp,  the 
P  station  residual 

•  -  -  T52*  Stp)  ' 

Figure  5  shows  o  as  a  function  of  6tp  for  6tp  ranging  from 
0  to  3  seconds,  for  a0=  8.2  km/sec  and  various  values  of  D. 

It  should  be  noted  that  the  starting  point  for  the  calculation 
corresponds  to  a  high  velocity  upper  mantle,  and  thus  to  the 
early  arrivals  rather  than  to  the  zero  P  residual  of  Cleary- 
Hales  (1966)  . 


THE  ELASTIC  CONSTANTS  AND  DENSITY 


If  low  shear  velocities  in  the  uppej  mantle  are  due 
to  an  approach  to  melting,  perhaps  only  partial,  then  the 
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principal  effect  would  be  cn  the  shear  modulus  y.  It 
was  suggested  (Hales,  1964)  that  in  such  cases  6a/a  would 
lie  between  2/3  and  4/9  of  68/6.  The  upper  limit  was  derived 
from  a  model  in  which  A  remained  constant,  and  the  lower  from 
one  in  which  k  remained  unchanged.  These  limits  correspond  to 
6 6/ 6a,  lying  between  0.87  and  1.30,  wheveas  the  observed 
value  is  1.24,  the  95%  confidence  limits  being  0.96  and  1.53. 
Thus  the  model  in  which  k  is  constant  fits  the  observations 
reasonably  well.  No  allowance  has  been  made  for  the  effect 
of  possible  changes  of  density. 

The  attenuation  of  seismic  energy  in  the  upper  mantle 
was  discussed  by  Anderson,  3en~Menahtm  and  Archambeau  (1965) . 
They  found  that  the  attenuation  of  surface  waves  could  be 
represented  by  a  model  in  which  the  elastic  constants  were 
complex,  the  imaginary  component  of  y,  y*  being  relatively 
large  in  the  low  velocity  zone.  They  concluded  that  a  model 
in  which  the  imaginary  part  of  the  incompressibility  k*  =  0, 
i.e.  in  which  there  were  no  losses  in  pure  compression, 
seemed  to  be  a  good  cne,  although  experimental  error  permitted 
values  of  k*  as  high  as  y*/2.  Thus  the  attenuation  studies 
suggest  also  that  the  changes  in  y  are  larger  than  those  in 
the  other  elastic  constants. 

There  is  good  evidence  that  there  are 
density  anomalies  in  the  upper  mantle  to  the  west  of  the 
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Rocky  Mountain  front.  Seismic  refraction  studies  by  the 
Crustal  Structure  Branch  of  the  U.  S.  Geological  Survey 
(Pakiser  and  Zietz,  1965)  show  chat  the  continental  crust  in 
the  western  United  States  is  generally  thinner  than  40  km, 
whereas  in  the  central  and  eastern  United  States  the  crustal 
t!  ickness  is  generally  greater  than  40  km.  Although  some 
structures  such  as  the  Sierra  Nevadas  arc  locally  compensated, 
it  is  clear  that  a  considerable  part  of  the  compensation 
»  for  the  elevated  region  to  the  west  of  the  Rocky  Mountain 

front  comes  from  lower  densities  in  the  upper  mantle. 

The  Bouguer  gravity  anomalies  in  the  western  plateau 
region  are  of  order  -200  milligals  over  large  areas  (Woollard 
and  Joesting,  .964)  .  It  is  not  clear  how  much  of  these 
large  negative  anomalies  should  be  ascribed  to  density 
variations  in  the  upper  mantle,  but  200  milligals  would  appear 
to  be  an  upper  limit. 

Supposing  that  the  densit /  anomaly  is  6p,  extends 

over  D^km  depth  range,  and  can  be  represented  as  an  infinite 

/. 

plate,  then  D6p  =  4.75  and  =  1.4  km.  For  arrivals  at 

p  o 

teleseismic  distances , the  rays  in  the  upper  mantle  are  inclined 
at  angles  of  about  30°  to  the  radius,  so  that  D'  =  0.9  D  approx¬ 
imately  an^  thus  D6p  is  5.2.  Since  the  range  of  S  station  anomaly 

—  =  P  sec,  is  equal  to  36  km.  Since  6  =  it  follows 

8  2  So  *  P 

o 

that  6  6/6  =  *2—  -  *5—,  and  thus  the  contribution  to  6  6/S0 

0  N  p0 

of  the  density  change  is  about  1/24  of  the  contribution  from 

6 u/u  .  Furthermore  it  is  of  opposite  sign. 

0 
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The  model  of  elastic  constant  change  which  best  fits 
the  observations  is  one  in  which  y  alone  changes  and  k 
remains  constant.  For  this  model  the  ratio  of  the  S  residual 
to  the  P  residual  is  3.9,  whereas  the  observed  value  was  3.72. 
The  95%  confidence  limits  of  the  observed  value  are  2.86  and 
4.58.  If  we  fit  a  model  to  the  l^>wer  limit  of  2.86,  dk/k0 
is  found  t.o  be  5y/4y  approximately. 

In  summary  then,  the  observed  P  and  S  travel  time 
anomalies  are  consistent  with  ?  model  in  which  y  alone 
changes,  k  remaining  constant,  but  relative  charges  in  k  up 
to  25%  of  those  in  y  are  allowed  by  the  observational  data. 

For  a  model  in  which  X  remains  constant  the  ratio  of  the 
S  anomalies  to  the  P  anomalies  is  2.60  and  lies  outside  the 
95%  confidence  limits  for  the  observations. 


Birch  (1961)  has  given  a  relation  between  a  and  p 

for  a  wide  range  of  rock  types.  From  Birch's  data  it  follows 

that  da/dp  is  about  3  km/sec  per  gm/cm3  for  most  rocks.  Since 

6tp  has  a  range  of  2’  seconds  and  is  equal  to  ,  it  follows 

that  D<5u  -  134.4  kmr/sec.  it  has  been  3hovn  that  DSp  is  of 

order  5.2.  Thus  -j—  =  26  and  is  an  order  of  magnitude  larger 

than  the  value  expected  if  Birch's  relation  were  applicable 

in  this  region  of  the  upper  mantle.  The  infinite  plate 

calculation  used  above  may  not  be  appropriate.  For  reasonable 

three  dimensional  structures  Dop  lies  between  10  and  15  and 

.12.  between  13  and  9,  considerablv  greater  than  value  of  about 

<5p 

3  expected  of  Birch's  relation. 
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THE  EFFECT  OF  TEMPERATURE 

Experimental  determinations  of  the  effect  of  temperature 
on  the  velocity  of  shear  waves  in  rocks  at  high  pressure 
were  given  by  Birch,  Schairer  and  Spicer  (1942,  p.  83)  and 
Birch  (1963).  For  ultrabasic  rocks  the  relative  change  in 
6  is  of  order  0.004  to  0.008  per  100°  C.  Since  the  range 

r  q 

of  the  S  residuals  is  8  seconds,  -5—  =  0.378,  0.189,  0.0756 

Bo 

and  0.038  for  D  =  100,  200,  500  and  1000  km  respectively. 

Thus  if  the  difference  in  temperature  were  spread  over 

1000  km  depth  range,  the  temperature  difference  required  would 

be  500°  to  1000°  C. 

However,  Anderson,  Ben-Menahem  and  Archambeau  (1965) 
found  low  values  f  Q  over  a  depth  range  of  about  300  to  400  km, 
so  that  if  the  reduction  in  shear  velocity  and  the  low  Q 
values  are  associated,  as  seems  probable,  the  major  part  of 
the  effect  would  occur  over  300  to  400  km,  and  considerably 
greater  temperature  differences  than  500°  C  would  be  required. 

Soga,  Schreiber  and  Anderson  (1966)  have  estimated  the 
seismic  velocities  at  high  temperatures  from  those  at 
relatively  low  temperatures.  For  Mg2SiO/j  9y/9T  is  about 
0.0060  per  100°  C.  Thus  using  the  Soga  et  al .  data  the 
temperature  changes  required  would  be  of  order  600°C  over 
1000  km  and  1250°C  over  500  km. 

The  Soga  et  al.  equations  are 
a  =  7.750  -  0.000362T  and 

3  =  4.513  -  0.000271T  to  the  first  order  for  Mg2Si04 
where  T  is  the  temperature.  From  these  equations  it  can  be 

r  q  2 

deduced  that  for  changes  of  temperature  is  2.25, 
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which  is  significantly  less  than  the  observed  value  of  6tg/6t 
It  is  therefore  unlikely  that  the  S  wave  station  anomalies 
can  be  accounted  for  in  terms  of  changes  in  temperature 
unless  the  temperatures  are  such  that  one  component  of  the 
system  approaches  its  melting  point  or  some  new  thermoelastic 
effect  occurs. 

If  one  component  of  the  system  were  to  approach  melting 
point,  or  if  the  intercrystalline  constraints  were  relaxed  in 
some  other  process,  changes  in  shear  velocity  much  greater 
than  found  by  Birch  and  Bancroft  (1938,  1940),  Birch  (1963) 
or  Sega  et  al.  (1966)  nr  ght  occur. 

CORRELATION  OF  STATION  RESIDUALS  WITH  HEAT  FLOW 

It  was  pointed  out  by  Cleary  and  Hales  (1966)  that 
station  anomalies  were  negative  in  shield  areas  and  that 
these  were,  in  general,  regions  of  low  heat  flow.  Recently 
published  heat  flow  values  by  Hyndman  (1967)  for  eastern 
Queensland  average  1.2  y  cal/cm2  sec,  whereas  those  in 
southeastern  Australia  and  Tasmania  have  a  mean  of  2.36 
±  0.38  y  cal/cm2  sec  (Howard  and  Sass,  1964).  The  P  station 
anomalies  in  southeastern  Australia  and  Tasmania  are  positive 
whereas  those  for  Brisbane  and  Charter' s  Towers  in  eastern 
Queensland  are  negative. 

Thus  there  appears  to  be  a  correlation  between  positive 
station  resi aals  and  high  heat  flow.  It  is  our  opinion 
that  the  station  anomalies  originate  largely  in  the  upper 
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mantle,  and  so  it  would  appear  that  the  high  heat  flows 
are  associated  with  higher  temperatures  in  the  upper  mantle 

CONCLUSION 

The  general  conclusion  ao  be  drawn  from  this  discussion 
of  P  and  S  travel  time  anomalies  in  the  United  States  is  that 
the  differences  between  the  central  and  western  regions  are 
due  in  the  main  to  higher  temperatures  in  the  upper  mantle 
of  the  western  United  States  such  that  there  is  an  approach 
of  one  constituent  to  melting  with  a  consequent  reduction 
of  a,  or  alternatively  that  the  intercrystalline  constraints 
are  somehow  reduced  at  high  temperature. 

There  is  no  doubt  that  some  parts  of  the  travel  time 
deviations  arise  from  differences  of  crustal  structure  and 
differences  of  composition  ir.  the  crust,  but  this  study 
suggests  strongly  that  the  major  part  of  the  deviations  for 
the  two  regions  are  to  be  ascribed  to  differences  in 
temperature  conditions  in  the  upper  mantle.  The  scatter  in 
the  ratio  6tg/6tp  is  probably  due  in  some  measure  to  the 
differences  in  crustal  structure  for  which  it  is  probable 


that 


■5? 
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would  be  of  order  1.7  to  1.8. 
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Figure  Captions 


Figure  1.  P  and  S  station  anomalies 

Figure  2.  a,  a/8  and  a  as  functions  of  8  for  a  model  in 
which  a0  =  8.2  km/sec,  8  0  =  «„//?. 

Figure  3.  Differences  between  the  calculated  travel  times 

for  various  models  and  those  for  the  Gutenberg 
velocity  distribution  (Model  0) .  The  dashed 
curve  gives  differences  between  calculated 
ScS  travel  times  for  the  McEvilly  Love  wave 
model  and  those  for  Model  0. 


Figure  4.  Differences  between  calculated  travel  times  for 
several  models  and  the  Model  0. 

Figure  5.  o  as  a  function  of  6tp  for  6tp  ranging  from  0  to 

3  sec,  for  ctc  =  8.2  km/sec  and  various  values  of 
D. 


Table  Captions 

Table  1.  '"ources  of  models. 

Table  2.  Velocity  distributions  for  the  models.  The  bars 
indicate  discontinuities. 
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TABLE  1.  SOURCES  OP  MODELS 


Model  No. 


4 

5 

6 
8 

10 


Gutenberg:  from  the  table  in  Handbook  of 

g^ysical  Constants  (1966)  compTied  by  P.  Press, 
standard  crust  added.  Times  for  this  model  were 
used  as  bases  for  Figures  3  and  4. 

Jeffreys:  from  the  table  in  Handbook  of  . 

Physical  Constants  (1966) , "compiled 'by  F  Press 
07ust  with  velocity  3.6  km/sec  added/  Press- 


Models  4,  5,  6,  8  and  10  have  Gutenberg  lower  mantle, 
and  increased  velocities  in  the  upper  mantle  as 
compared  with  model  0. 


2 

11 

3 


12 


Anderson  model  (ghysics  and  Chemistry  of  the  Earth, 
at  1400^  196°'  ^ted  to  Gutenberg  lower  mantle 

McEvilly  (1964)  Rayleigh  model  fitted  to  Gutenberg 
lower  mantle  at  500  km.  9 

u°Ve  model  fitted  to  Gutenberg  lower 
mantle  at  500  km. 

Alexander  (1963)  crust  and  upper  mantle  fitted  to 
Gutenberg  lower  mantle  at  500  km. 

Brune  and  Dorman  (1963)  Canadian  Shield  model  plotted 
to  Gutenberg  lower  mantle  at  500  km. 

MacDonald  and  Ness  (1961)  Gutenberg  IV  model.  Note 

Asa  result0?^  S?~iStS  °f  iayers  °f  constant  velocity. 
As  a  result  the  differences  from  the  Gutenberg  Handbook 
model  are  irregular. 
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